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Abstract 



Max cones are max-algebraic analogs of convex cones. In the present 
paper we develop a theory of generating sets and extremals of max cones 
r^ ' in R" . This theory is based on the observation that extremals are minimal 

elements of max cones under suitable scalings of vectors. We give new 
proofs of existing results suitably generalizing, restating and refining them. 
Of these, it is important that any set of generators may be partitioned 
into the set of extremals and the set of redundant elements. We include 

rsl ' results on properties of open and closed cones, on properties of totally 

^ ' dependent sets and on computational bounds for the problem of finding 

^^ , the (essentially unique) basis of a finitely generated cone. 



in 

o 

^ ■ 1 Introduction 
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j^ , By max algebra we understand the analog of linear algebra obtained by consid- 

ering R+ (the nonnegative reals) with max times operations: 



a (Bb := max(a, 6), 

a (^ b := ab 

extended to matrices and vectors. That is, if A = {aij), B = (bij) and C = (cij) 
are matrices of compatible sizes with entries from K+ and a G M+, we write C = 
AQ)B iicij = aij®bij for all i, j, C = A(8)i? if Cy = 0j, aik®bkj = ma.Xk{aikbkj) 
for all i , j and C = a (i) A ii Cij = a®aij for all i,j. There are several essentially 
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equivalent ^ definitions of max algebra. An early paper presenting the above 
version is |Vor| . another early paper presenting an equivalent version is |CGO| . 
For more information on max algebra, its generalizations and applications the 
reader is referred e.g. to |B(^0Q| , [Eut], ^nnT\ \CG2 CGS], [HMl and [UZ]- See 
also |LM| for recent developments in the area and for further references. 

We give a summary of the contents of this paper. In Section 2 we begin by 
defining (max) cones, extremals, generating sets, independent sets and totally 
dependent sets. The key observation is Proposition ^] which extends [Josl 
Proposition 2.9]. We deduce in Theoremll4lthat extremals are minimal elements 
of max cones under suitable scalings of vectors. This leads us to a reformulation 
and new proof of the basic result Theorem 1161 which is also easily derived from 
|Wagl| Proposition 2.5.3]: Every generating set of a max cone can be partitioned 
into the set of the extremals of a cone and a set of redundant elements. It 
follows that if a cone has a basis then (under a scaling condition) it consists of 
the extremals of the cone and hence the cone has a basis unique up to scalar 
multiples, see Theorem ^| and its corollaries for more detail. In Corollary 
|55] we then turn to the case when the set of extremals of the cone is empty, 
in which case every generating set is totally dependent. Discussion of totally 
dependent sets specifically may be new. Towards the end of this section we 
consider topological notions. In Corollary 1231 we show that (under a natural 
restriction) every open cone has totally dependent generators and in Proposition 
1241 we prove an analogue of Minkowski's theorem for closed cones. This result 
extends a result due to |Hel| and it also appears as |GK1I Theorem 3.1] where 
a different proof is given. 

In Section 3 we give two simple versions of an algorithm, based on JCG2I The- 
orem 16.2] for finding the (essentially unique) basis of a finitely generated max 
cone and a MATLAB program which implements one version. We also relate 
our problem to the classical problem of finding maxima of a set of vectors which 
is described in |KLP | and in jPSL Section 3], and give the bounds for computa- 
tional complexity. 

We now relate the concepts and techniques of our proofs to those in other pub- 
lications. Most of our concepts appear in |Hel| . |WagO| and |Wagl| , sometimes 
under different names. For instance, extremals are called irreducible elements 
in |WagO and |Wagl| , and minimal elements are called efficient points in |IIel| . 
Our key Proposition II II mav also be derived from (possibly slightly extended) 
results found in some of our references. Examples are results in |Vor| Section 2 
in terms of set coverings, see also |CGO| . |CG2I Theorem 15.6] and |Butj Section 
2, or the fundamental results of |CG2| and |CG3[ Section 3] concerning max lin- 
ear systems. The latter are also found in [CGlj and |US| in terms of a projection 
operator. Further, such results on max linear systems as Proposition 1111 and 
Proposition!^ can be extended to the case of functional Galois connections, as 
it is shown in |AGK| . The generalizations considered in | AGK| are useful in 

^that is, algebraically isomorphic 



many areas including abstract convex analysis, the theory of Hamilton- Jacobi 
equation and the Monge-Kantorovitch mass transportation problem. 

Our topic is also related to (and partially stimulated by) the emerging field of 
tropical geometry which develops basic concepts of max algebra in a different 
form and with different terminology and applies these to finitely generated struc- 
tures, see |DS | .| Jos) and |B Yuj . In particular, Proposition^^can also be seen as 
a minor extension of |Josl Proposition 2.9] which is important in the theory of 
tropical halfspaces. The emphasis of these papers is on geometry, while in this 
paper it is on algebraic and order theoretic results. Max cones are also studied 

in | GKO| and |GK1| . The main effort of these papers is to develop the theory 
of max-plus convex sets and their recession cones. This theory is not present 
in our paper. In turn, we deal with more general cones and we emphasize the 
link to set maxima and give a more detailed description of bases and generating 
sets. 

In max algebra as in linear algebra a basis is normally defined as an indexed set, 
that is a sequence if the basis is finite or countable, see |Wagl| for a definition 
in max algebra or |Boul p. 10] in linear algebra. Since we wish to show the 
inclusion of the set of extremals (which do not have a natural order) in every 
generating set or basis for a cone we define the latter in term of sets in Section 
2. We thereby exclude the possibility of a repetition of elements in generating 
sets. But we change our point of view in Section 3 on algorithms for finitely 
generated cones since we wish to consider the generators as columns of a matrix. 



Max cones have much in common with convex cones, see |Roc| for a general 
reference. This has been exploited (and generalized) in many papers including 
those just quoted and e.g. |CGQS| and KZ . To this end, the basic concepts of 
this paper and such results as Theorem 1161 Prooosition 1241 and Proposition 1251 
have their direct analogs in terms of positive linear combinations and in convex 
analysis. We do not provide details, as convex geometry is also beyond our 
scope here. 

2 Generating sets, bases and extremals 

We begin with two standard definitions of max algebra. 

Definition 1 A subset K of R" is a max cone in R" if it is closed under © and 
(g) by nonnegative reals. 

Definition 2 Let S C M" . Then u is a max combination of S if 



u = 

xes 



0A,a;, A, eR+, (1) 



where only finite number of Aj, ^ 0. The set of all max combinations will be 
denoted by span(S'). We put span(0) = {0}. 



Evidently, span(S') is a cone. If span(S') = K, we call S a set of generators for 
K. 

Definition 3 An element u € K is an extremal in K if 

u = V (B w, v,w ^ K =^ u = V ov u ~ w. (2) 

If u is an extremal in K and A > then Am is also an extremal in K . 

Definition 4 An element x G M" is scaled if ||a::|| — 1. 

For most of this section, ||a;|| may be any norm in R" (they are all equivalent). 
However, in the end we specialize to the max norm, ||a;l| — maxxi, in order to 
exploit the property that it is max linear on R". If 5 C R" we may call S 
scaled to indicate that it consists of scaled elements. 

Definition 5 Let 5* be a set of vectors in R" . 

1. The set S is dependent if, for some x G 5, a; is a max combination of 
S \ {x}. Otherwise, S is independent. 

2. The set S is totally dependent if every a; G S* is a max combination of 
S\{x}. 

Thus the empty set of vectors is both independent and totally dependent. Since 
span(0) = {0}, the set {0} is totally dependent. 

Definition 6 Let K he a, cone in R" . A set S of vectors in R" is a basis for K 
if it is an independent set of generators for K. 

The set of all unit vectors {e^, p = 1, . . . , n} defined, as usual, by 



e?= ; /' (3) 




is a basis of R" , which is called standard. 

Lemma 7 Let S he a set of scaled generators for the cone K in R" and let u 
he a scaled extremal in K . Then u ^ S. 

Proof. Suppose u is given by the max combination (^). Since the number of 
nonzero A^; is finite, we may use DefinitionOand induction to show that u = XxX 
for some x. But u and x are both scaled, hence u = x and u € S. U 

Lemma 8 The set of scaled extremals of a cone is independent. 

Proof. If the set E of scaled extremals is nonempty let w be a scaled extremal 
in K and apply Lemma[71to the cone Ki :— spa,n{E\{u}) . This shows u ^ Ki 
and the result is proved. D 

Below we use subscripts for elements of vectors in R" and superscripts to 
label vectors. 



Definition 9 Let v G R" . Then the support of v is defined by 

supp(t;) = {j e {1, . . . , n} : Vj > 0}. 
The cardinahty of supp(u) wiU be written as | supp(u)|. 

In order to relate the natural partial order on R" to results on extremals of 
cones we introduce a scaling of vectors in M" for each j G {1, . . . , n} such that 
for each scaled vector Vj — 1. 

Definition 10 1. Let u G M" and suppose j G supp(u). Then we define 

u{j) = U/Uj. 

2. Let S C El:. We define S{j) = {u{j) : u G S" and j G supp(m)} for all 
j = l,...,n. 

3. Let S C M" . An element u G S" is called minimal in S, ii v < u and v € S 
implies that v — u. 

4. Let K he a, cone in K" , let u G K, and let j G supp(w). We define 

Dj[u) = {veK{j):v<u[j)}. 

Our key observation is the following Proposition. It can be viewed as a minor 
but needed extension of |Jos[ Proposition 2.9], see also the remarks concerning 
it in our Introduction. 

Key Proposition 11 Let S C M" . Then the following are equivalent: 

1. u ^ span(S'). 

2. For each j G supp(u) there is an x^ G S such that j G supp(a;-') and 

x^ij)eD,iu). 

Proof. 2. =^ 1: If 2. holds, then u = ®j^supp(u) ^J^'' where A^ — Uj/x-j. 
1. =^ 2. Conversely if 1. holds, then it follows immediately from QJ that for 
each j G supp(u) there is an x^ G 5 with XjX^ < u and {XjX^)j = Uj. Clearly, 
Xj — Uj/x-j which yields 2. D 

The following immediate corollary to Proposition II II is essentially found as 
jHell Theorem II. 1] and as jDSI Proposition 5]. It is analogous to Caratheodory's 
Theorem. 

Corollary 12 Let S C R" . Then u G span(iS') if and only if there are k vectors 
x^, . . . , x*^ G 5, where k < \supp{u)\, such that u G spanja:^, ...,x^} . D 

Corollary 13 Let K he a cone in K" and let T be a set of generators for K . 
Let U Q T and let S = T\U . Then S generates K if and only if each u ^ T 
satisfies condition 2. of Proposition MIX D 



Theorem 14 Let K be a cone in M" generated by S and let u £ S, u ^ 0. 
Then the following are equivalent: 

1. u is an extremal in K. 

2. For some j G supp(u), u{j) is minimal in K{j). 

3. For some j G supp(u), u(j) is minimal in S{j). 

Proof. 1. =^ 3. If I supp('u)| = 1 then u{j) is minimal in S{j). So suppose that 
I supp(u)| > 1 and that u(j) is not minimal in S{j) for any j G supp(m). Then 
for each j G supp(u) there exists x^ G S{j) such that x^ < u{i), x^ ^ u[j). 
Therefore u = ®,gsupp(M) '^j^'' ^ ^-nd u is proportional with none of x^ . Hence u 
is not an extremal in K. 

3. =^ 2. Let V € K and assume that j G supp(i;) and v{j) < u{j). We need to 
show that v{j) — u{j). By Proposition ^] there is a w G S" such that w{j) < 
v{j). Thus w{j) < v{j) < u{j) and by 3. it follows that w{j) = v{j) = u{j). 
2. =^ 1. Let u{j) be minimal in K{j) for some j G supp(u) and suppose that 
u = V (B w, v,w E K. Then both v < u and w < u and either Vj — Uj or 
Wj — Uj, say (without loss of generality) that Vj = Uj. Hence v{j) < u{j) and 
it follows from 2. that v{j) = u{j). Hence also v = u which proves 1. 
D 

Note that in Theorem E| we can of course have S — K. Also note that 
CoroUarv 1131 mav be combined with Theorem II 41 to yield conditions for a set of 
generators to be redundant. 

Corollary 15 Let K be a cone in R" . If Dj{u) has a minimal element for each 
u E K and each j G supp(u), then K is generated by its extremals. 

Proof. Suppose that x^ is a minimal element oi Dj(u). Since, for v G K{j), v < 
x^ implies that v G Dj{u), x^ is also a minimal element of K{j). We now obtain 
the Corollary by combining Proposition II II and Theorem ^1 D 

Essentially, the following fundamental result was proved in |Wagl| Proposi- 
tion 2.5.3]. We suitably restate it: every set of generators S for a cone K can 
be partitioned as -E U i^, where ii^ is a set of extremals for K and the remainder 
F is redundant. Our proof is a combination of Proposition llll and Theorem ll4l 

Theorem 16 Let S be a set of scaled generators for a cone K in R" and let 
E be the set of scaled extremals in K. Then 

1. ECS. 

2. Let F — S\E. Then for any u E F , the set S\{u\ is a set of generators 
for K. 



Proof. Assertion n repeats Lemma [7| 

To prove Assertion [2 let u G F. Since u is not an extremal, by Theorem 1 141 for 

each j G supp(u) there is z^ G K such that z^ [j) < u{j). Since K = span(S'), by 
Proposition ^2 we also have y^ e S such that y^{j) < z^{j) < u{j). Evidently 
y^ ^ u, and applying Proposition !! II again, we get that u is a max combination 
of {y^ : j G supp(u)}, where y^ € S are different from u. Thus in any max 
combination involving u, this vector can be replaced by a max combination of 
vectors in S'\{u}, and the theorem is proved. D 

The following example shows that the set F of Theorem ^] need not be 
totally dependent. 

Example 17 Let K be the cone in i?^ generated by u^ = [1, l/r]^ , r = 1, . . .. 
The elements of K scaled with respect to the max norm are [1, a]'^ with < a < 

1. Thus u^ is the unique scaled extremal in K. But the set F — {u^ : r = 2, . . .} 
is not totally dependent since u^ is an extremal in span(i^) whose scaled elements 
are [1, a]^ with < a < 1/2. 

The following is a refinement of Theorem^] and also of |Wagl| Theorem 5]. 

Theorem 18 Let E be the set of scaled extremals in a max cone K . Let S G K 
consist of scaled elements. Then the following are equivalent: 

1. The set S is a minimal set of generators for K. 

2. S = E and S generates K . 

3. The set S is a basis for K . 

Proof. 1. =^ 2. By Theorem 1161 we have S = E U F where every element of 
F is redundant in S. But since 5* is a minimal set of generators, we must have 
^ = 0. Hence S ^ E. 

2. =^ 3. The set E is independent and a generating set. 

3. =^ 1. By independence of S the span of a proper subset of S is strictly 
contained in span(S'). D 

Theorem ^| shows that if a cone has a (scaled) basis then it must be its 
set of (scaled) extremals, hence the basis is essentially unique. We note that a 
maximal independent set in a cone K need not be a basis for K as is shown by 
the following example. 

Example 19 Let K C M^ consist of all [xi, X2Y' with xi > 2:2 > 0. If 
1 > a > 6 > 0, then {[1, a]-^, [1, b]^} is a maximal independent set in K which 
does not generate K. 

The following corollary is found e.g. as [DHI, Proposition 21], 'Jos', Proposition 
2.5] and also in Ser. Proposition 1], where it is used to obtain uniqueness results 
for definite max-plus matrices. As a special case of this corollary, the standard 
basis of W] is essentially unique. 



Corollary 20 If K is a finitely generated cone, then its set of scaled extremals 
is the unique scaled basis for K . 

Proof. Since K is finitely generated, there exists a minimal set of generators S. 
By Theorem [THIS' = E and S" is a basis. D 

Note that in the tropical geometry |DS| and |Jos | vertices of a polytope 
are defined to be the essentially unique generators determined in Corollary I2UI 
and hence vertices correspond to our extremals (and to Wagneur's irreducible 
elements). Next we obtain some corollaries concerning totally dependent sets. 

Corollary 21 If S is a nonempty scaled totally dependent set in M!f. then S is 
infinite. 

Proof. Suppose that S is finite and let K = span(S'). By Corollary 1201 K 
contains scaled extremals which, by Theorem 1 161 must be contained in S given 
that K = span(S'). But then S is not totally dependent. This contradiction 
proves the result. D 

Corollary 22 Let K he a cone in M" . The following are equivalent: 

1. There is no extremal in K . 

2. There exists a totally dependent set of generators for K. 

3. Every set of generators for K is totally dependent. 

Proof. Since there always exists a set of generators for K (e.g. K itself), each 
of the Conditions 2. and 3. is equivalent to Condition 1. by Theorem 1161 D 

We now consider W^ in the topology induced by the Euclidean topology of 
R". That is, a set in R" will be called open if and only if it is the intersection 
of an open subset of R" with R" . A cone K is called open if iir\{0} is open, 
and it is called closed if it is closed as a subset of R" , or equivalently of R". 

Corollary 23 // K is an open cone in R" that does not contain unit vectors, 
then every generating set for K is totally dependent. 

Proof. It is enough to show that there is no extremal in K , for then the result 
follows by Theorem^l Let u ^ K. Since u is not a unit vector, there are at least 
two indices fc, I £ supp(w). Since K is open, we have w^ = u — ee^ ^ K, p = k,l 
for sufficiently small e and u — w'' (B u/ . None of w^ , p = k,l is equal to u, 
hence u is not an extremal, and the corollary follows. D 

An example of an open cone in R" is furnished by the cone K of all positive 
vectors in R". We note that, for this particular case, CoroUarv 1231 was shown 
in |UGB| . Another example of an open cone consists of all vectors [a, bf" 
in R^ with a > 6 > 0. We acknowledge the following Proposition, which is 

analogous to Minkowski's Theorem, to [GKll Theorem 3.1]. There the result is 
proved directly by a minimality argument; here we deduce it from a corollary 
to Theorem 1141 which characterizes extremals of cones that may not be closed. 
It extends earlier results of IHell . 



Proposition 24 Let K be a closed cone in M". Then K is generated by its 
set of extremals, and any point in K is a max combination of not more than n 
extremals. 

Proof. Let u G K and let j G supp(u). It is easily shown that Dj{u) is compact 
since K is closed. Hence Dj{u) contains a minimal element x^ . The result now 
follows by CoroUarvllSland Corollary 1 121 D 

The max norm is max linear on M" : 

||Au©^t;|| = A||u||©/i||i;||. (4) 

This is exploited in the following proposition. 

Proposition 25 If S C R!J. is compact and ^ S, then the cone K = span(5') 
is closed. 

Proof. Consider a sequence u^ G K converging to v. Then, by Corollary 1121 we 
have 

n 

where w'" G S and A,^ £ M+. By gj) 

IKIH0A..IK1I. (5) 

Since the sequence u* converges (to v), the norms ||m*|| are bounded from above 
by some Mi > 0. On the other hand, we have ||w**|| > M2 for some M2 > 0, 
since S is closed and does not contain 0. Then by (|3J| Ais||?i;**|| < Mi for all i 
and s, and Xis < AI1M2 for all i and s. Thus A^s are bounded from above. 
But ||w*^|| are also bounded from above, since S is compact. This implies that 
there is a subsequence w^^'^ such that for all ,s = 1, . . . , n the sequences w-'(')* 
and Xjiijs converge. Denote their limits by w* and As, respectively, then w^ G S* 
and As < MiMj"^. By continuity of © and (8> we obtain that 

n 

^;-0AsU;^ (6) 

Thus V e K. n 

Corollary 26 // the set of scaled extremals of a max cone K is closed and 
generates K , then K is closed. 

Corollary 27 Any finitely generated max cone K is closed. 

We now give a counterexample to the converses of Corollary E] and Proposi- 
tion!^ (part 1), and to the converse of Corollary [5^1 (part 2). 



Example 28 1. In M.\ let S consist of all vectors [xi, X2, 1]^, < xi < 1/2 
such that a;i + X2 = 1 and let K = span(5'). Then the section of K given 
by 2^3 = 1 consists of all vectors [xi, X2, 1]^, < xi < 1/2, < X2 < 1 
such that xi + a;2 > 1. Note that S is the set of extremals of K scaled 
with respect to the max norm, but K = span(5') is not closed and for any 
u E K there are no minimal elements in Di(u) and D2{u). 

2. Now let S" = 5 U {u}, where u = [1/2, 0, 1]^ and let K' ^ span(5'). 
Then the section of K' given by 0:3 = 1 consists of K together with the 
line segment whose end points are u and [1/2, 1, 1]-^. Thus K' is closed. 
The set of scaled extremals of K' is 5" which is not closed. 

The cross sections of K and K' by X3 = 1 are shown on Figure 1, together with 
the generating sets S and S" = S* U {u}. 
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Figure 1 : Max cones of Example EHl 



3 Algorithmic considerations 



As explained in the introduction we redefine our basic concepts for this section 
which is concerned with finitely generated cones. We also restate a suitable 
adaptation of CorollarvEHl 



Definition 29 Let V 
deleting column i, i = 



1, 



V]^ and let V^ be the matrix obtained from V by 
. . , fc. Then the cone K generated by the columns 



,v oi V consists of all vectors of form V ® x, a; G 



1%. 



Further, the 



columns of V form a basis for K if, for i, i 
that V~,®x = v\ 



1, . . . , fc, there is no a; G M^ ^ such 



Proposition 30 Let V £ 



Then there exists 



ibmatrix U G 



jnp 



,0 < 



p < k whose columns form a basis for the cone generated by the columns of V 
(and every other basis is of form UPD, where P is a permutation matrix and 
D is a diagonal matrix with nonzero diagonal elements). D 



10 



We shall apply the following proposition. Note that all statements in this propo- 
sition have been proved in a more general setting in |CG1| . See also |But| , |CG3[ 
Chapter III] and |Vor . 

Proposition 31 Let U G M"''' with all columns nonzero and let v £ MJ^. Let 
X G M'ji he defined by 

Xi — min{t)j/u* : u* / 0, j = 1, . . . , n} (7) 

for i — 1, . . . ,k. Then 



U®x < V, (8) 

X = max{z G M+ : t/ ® z < u}, (9) 

U®x = max{C/® z : z G IR+,C/(g)z < w}. (10) 



Further, there exists z G M*^ such that U ® z = v if and only if U ® x = v. 

Proof. Assertion ||HJ) follows from the observation that U ® z < v \i and only if 
Zi < Vj/u^j if j G supp(u'), i — 1, . . . ,k. Note that x G M^ since no column of 
U is zero. Since (8 is isotone (that is, x < y implies A (g) a; < A^ y), assertions 
(O and 11U|) follow immediately. For the final statement assume that U (E) z = v 
for some z. By ^ and ((Tn|l we have v = UiS)z<U(^x<v, and the statement 
follows. The converse is trivial. D 

Algorithm 32 Input: V G M"''. 

Output: An n x p submatrix U of V whose columns form the essentially unique 

basis for the cone generated by the columns of V. 

Step 1. Initialize U — V . 

Step 2. For each j = 1 , . . . , A: if u^ 7^ set v = u^ , and for each i ^ j compute 

Xi by (TJ, if u* 7^ 0, and set Xi = otherwise. If U~x = v, set u^ = 0. 

Step 3. Delete the zero columns of U . The remaining columns of U are the 

basis we seek. 

Remark 33 The restriction in PropositionlH II that each column U G M"'^ must 
have a positive element was imposed to avoid definitions for a/0, a > 0, or 0/0. 
The restriction is inessential in the sense that for general U G W^ we may define 
Xi by whenever u' 7^ and choose Xi arbitrarily in K^ whenever u' = 0. 
Then all assertions of the Proposition still hold, with exception of ©• It is 
possible to extend M" by adding a maximal element cxo so that still holds. 

We omit details and present the MATLAB program maxbas that implements 
Algorithm l32l but employs such an extension. We also give an example with some 
elements equal to 0. Note that in |CG2I Theorem 16.2] a related algorithm called 
A-test has been presented. It enables us to identify columns that are dependent 
on other columns of an n x /c matrix in 0{nk^) time. However, there is no 
discussion of bases in connection with this method in |CG2| . 



11 



7othe unique max times basis for the max col space of A 
7ofunction [B,f] = maxbas(A), 

7oB = the unique max times basis for the max col space of A 
7,f = indices of columns of B in A 
7oCalls maxpr, max multiplication of matrices 
function [B,f] = maxbas(A) 
[m,n] = size(A); B = A; t = max(max(A)); 
for j = l:n 

V = compl(j ,n) ; 

c = B(:,j); BB = B(:,v); warning ( 'of f ') , 
e = ones(l,n-l); C = c*e; x = min(C./BB)'; 
z= maxpr (BB,x) ; 

if abs(c-z) < t*eps, B(:,j) = 0; end, 
end u = max(B) ; f = find(u >t*eps) ; B = B(:,f); 
A = 

Columns 1 through 5 
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We note that a second form of the algorithm may be based on set covering 
condition 1)1111 below, which appears in |CG2I Theorem 15.6]. It can also be 
found in |Vor| Section 2 and in |But| Section 2 but only in the case when all 
vectors are positive. See also jAGKI Theorem 3.5] for an interesting functional 
generalization of this condition (and more). With v and U as in Step 2, denote 
by Ni the set {j : v{j) > u^{j)}- By Proposition^ v G span(u^, . . . , u™) if and 
only if 

m 

U N, = supp(«). (11) 

With X given by Q, we note that 

^^ f{j e supp(u') : Wj/uJ ^xj ifx, ^0, 
[0, ii Xi = 0. 

Thus Step 2 in Algorithm 1321 mav be replaced by 
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Step 2': For each j — 1, . . . ,k such that u^ ^ 0: set v — u^ and for each i ^ j 
compute Ni = {j: v{j) > u^{j)} according to ((T^ . if w' j^ 0, and set Ni — % 
otherwise. If Ui=^i ^i — supp(w), set u^ = 0. 

The version with Step 2' is also weU- known. It is implemented in the max-plus 
toolbox of Scilab, a freely distributed software. See |GSI Sect.III-B] for the 
documentation. 

Our algorithms are of complexity 0{nk^). 

If 5* is the set of columns of the matrix [/, then it follows from Theorem E] 
that a basis for the cone generated by S consists of the union of the n sets 
-^(i); j = 1, . . . ,n, where M{j) consist of the vectors minimal in S{j). The 
problem of finding all maxima (or minima) of k vectors in R" is considered 
in |KLP| . and also in jPSL Section 4.1.3], where it is dubbed the problem of 
Erehwon Kings. The computational complexity of methods developed in |KLP| 
and | PS| is bounded from above by 0(n^fc(log2 fc)"~^) + 0(fclog2 k), n > 2, see 
| PS[ Theorem 4.9] and [KLPL Theorem 5.2]. ^ To solve our problem we can 
apply these methods to each S{j),j — 1, . . . ,n separately. Taking into account 
that for each j we need 0{nk) operations to find the coordinates of essentially 
(n — l)-dimensional vectors in S{j), this yields an alternative method with 
complexity not smaller than 0{n^k) and not greater than 0{n^k{log2 k)"~^) + 
0{k log2 fc) , n > 3. This method may be preferred if logj k is substantially larger 
than n. 
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